THE CRITICAL BEHAVIOUR OF THE 
SPIN-3/2 BLUME-CAPEL MODEL IN TWO 

DIMENSIONS 



J. C. Xavier , F. C. Alcaraz 
Departamento de Ffsica 
Universidade Federal de Sao Carlos 
13565-905, Sao Carlos, SP, Brazil 

D. Pena Lara * and J. A. Plascak 
Departamento de Ffsica 
Universidade Federal de Minas Gerais 
30161-970, Belo Horizonte, MG, Brazil 

PACS numbers : 64.60Kw, 64.60Cn, 64.60Fr 



Abstract 

The phase diagram of the spin-3/2 Blume-Capel model in two di- 
mensions is explored by conventional finite-size scaling, conformal in- 
variance and Monte Carlo simulations. The model in its r-continuum 
Hamiltonian version is also considered and compared with others spin- 
3/2 quantum chains. Our results indicate that differently from the 
standard spin-1 Blume-Capel model there is no multicritical point 
along the order-disorder transition line. This is in qualitative agree- 
ment with mean field prediction but in disagreement with previous 
approximate renormalization group calculations. 
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1 Introduction 



The Spin-S Blume-Capel Model is a generalization of the standard Ising 
model with dynamics described by the Hamiltonian 

H=-JY. SiS j + Dj2sl (1) 

<i,j> i 

where the first sum runs over all nearest neighbor and the spin variables 
assume values —S, —5 + 1, 5, associated with each site i of the lattice. 
In Eq.(|I]) J is the exchange coupling and D is a single spin anisotropy pa- 
rameter. In the case where 5 = 1 this Hamiltonian was proposed originally 
by Blume and Capel |1J] for treating magnetic systems and has also been 
used in describing 3 He- 4 He mixtures 0. This 5 = 1 model was studied by 
a variety of methods such as mean field [|TJ, two-spin cluster [0, variational 
methods f|, constant coupling approximation ||, Monte Carlo simulations 
[P, [jj, finite-size-scaling on its transfer matrix || or quantum Hamiltonian 



[[]. |10| and renormalization group methods [JTTJ, [12|. It is well established that 
for dimension d > 2 the 5 = 1 model presents a phase diagram with ordered 
ferromagnetic and disordered paramagnetic phases separated by a transition 
line which changes from a second-order character (Ising type) to a first-order 
one at a tricritical point. More specifically in two dimensions the machinery 
coming from conformal invariance [0, [TJj] indicates that at this tricritical 
point the long-range fluctuations are governed by a conformal field theory 
with central charge c = 7/10 ||, 1C]. In this case all the critical exponents 



and the whole operator content of the model are obtained |I(J . 

For values of spin 5 > 1 however the situation is quite unclear with fewer 
results and contradictions among them. The mean field calculation (see 



also Iflql ) predicts different phase diagrams for integer or half-odd-integer 
spins. For spin 5 = 3/2 it gives a second order phase transition with no 
tricritical point and a separated first-order transition line which ends up in 
an isolated multicritical point. In contradiction with these results a calcu- 
lation [O based on a renormalization group introduced in [T7| give us for 



the two dimensional model a unique first-order transition line at low tem- 
perature which ends up in the second-order transition line at a tetracritical 
point. Similar results are also obtained by using a Migdal-Kadanoff real 
space renormalization group approach |18[ . 



Motivated by these contradictions we decided to study extensively the 
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two-dimensional S = 3/2 model by using conformal invariance and finite-size 
scaling, since these methods give us in the S = 1 model the most precise and 
conclusive results. In order to supplement our conclusions, in the S = 3/2 
model, we also perform Monte Carlo simulations in the region of the phase 
diagram where multicritical behavior might occur. 

The layout of this paper is as follows. In section 2 the transfer matrix of 
the model and the relations used in our finite-size studies are presented. In 
section 3 we discuss initially our results for the S = 1 model. Although the 
phase diagram is well established we believe that our numerical estimate for 
the tricritical point is the most precise in the present literature. In section 4 
our results for the spin 3/2 is presented and discussed. In section 5 we study 
the spin-3/2 model in its r-continuum formulation and compare its quantum 
Hamiltonian with another spin-3/2 quantum chain which is known to present 
multicritical point along the order- disorder phase transition line. In section 6 
we close our paper with a summary and conclusions of our main results, and 
finally in appendix A the details and methods of our Monte Carlo simulations 
for the spin-3/2 model are explained. 



2 Transfer Matrix and Conformal Invariance 
Relations 

The row-to-row transfer matrix T of the Hamiltonian ([]]) in a square lattice, 
with horizontal width N has dimension (2S + 1) N x (2S+1) N . Its coefficients 
< s[, s' N \T\sx, sjv > are the Boltzmann weights generated by the spin 
configurations {s\, ...,sjv} and {s[, s' N } of adjacent rows. In general, if 
we consider (|1]) with periodic or free boundary conditions in the horizontal 
direction, the matrix T will have all its elements non-zero, being a dense 
matrix. This difficulty restricts the numerical study for very small lattices. 
A possible way to reach larger lattices, with reasonable computational effort, 



is to consider the model in its r-continuum limit [T9] . In this case after a high 
anisotropic limit of its couplings in the vertical and horizontal directions we 
replace the transfer matrix by a simpler effective quantum Hamiltonian which 
is a sparse matrix. Most of the finite-size scaling studies exploring conformal 
invariance, for general models were done in this formulation. However in 
such approach the calculation of the transition lines, in terms of the original 
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coupling (isotropic) is not possible. Another possibility, which is suitable for 
our purposes, is to impose helical boundary conditions along the horizontal 
direction. In this case the transfer matrix can be written as 

f = t N , (2) 

where t is the transfer submatrix with elements given by 

JV-l 

< S h S^Vl; -■> S N >= ex P (^ ls l( s 2 + s' N ) - dt^sf) Y[ $y+i.«j> ( 3 ) 

i=l 

with t = ksT/J (ks is the Boltzmann constant) and d = D/J. We also 
see from ([|) and (Q) that although T is a dense matrix t is a sparse one, 
having only a fraction / = exp (—(N — 1) ln(25* + 1)) of non-zero elements. 
The Hamiltonian ([!]), like the standard Ising model has a global Z{2) sym- 
metry. The helical boundary condition does not break this symmetry, which 
translates in the commutation of t with the parity operator 

N 

i=i 

where Ri = I<g)l(g)---I®i2®l---®1, and R, located at Ith position in the 
product, are (2S + 1) x (2S + 1) matrices given by 

s 

R= ^2 \s >< —s\. 

s=-S 

Consequently in the basis where P is diagonal we can separate the vector 
space associated to T into two disjoint sectors labelled by the eigenvalues 
p — ±1 of the parity operator. 

Let us for convenience introduce the Hamiltonian matrix H = — ln(T) = 
—Nhx(t), and denote by the e/ jP the Ith lowest (/ = 1,2,...) eigenvalue in 
the sector with parity p. The eigenenergies of the ground state and first 
excited state are given by e^i and e\-i, respectively. These two leading 
energies, and consequently the mass gap Gj\[(t, (i) = e x i - ex a, for a system 
of width iV are real and can be obtained for relative large iV by a direct 



application of the power method |20]. The standard finite-size scaling 



gives an estimate of the phase transition curve t c = t c (d). This curve is 



4 



evaluated by extrapolations to the bulk limit (N — > oo) of sequences t c (d, N) 
obtained by solving 

G N (t c )N = G N+1 (t c )(N + l) N = 2,3, ... . (4) 

An heuristic method, which was proved to be effective in obtaining multi- 
critical points in earlier works [§, |22|, is to simultaneously solve for 
three different lattice sizes 

G N (t c )N = G N+1 (t c )(N + 1 ) = G N+2 (t c )(N + 2) N = 2, 3, ... . (5) 

Once the transition curve is estimated, in the region of continuous phase 
transition we expect the model to be conformal invariant. This symmetry 
allows us to infer the critical properties from the finite-size corrections of the 



eigenspectrum at t c [13, 14]. The appropriate relations which give the scaling 



dimensions and conformal anomaly in the case of periodic or free boundary 



conditions are well known in the literature [13. 14]. In the present case, where 
we impose helical boundary conditions, these relations are slightly different 
from the periodic case since translational invariance is lost in the horizontal 
direction. The conformal anomaly c, which labels the universality class can 
be calculated from the large- A?" behavior of the ground-state energy 

-^ = ^- — + 0^ ), (6) 

where eoo is the ground-state energy, per site, in the bulk limit and v s , which 
is unity in our transfer matrix calculations, is the sound velocity. The scaling 
dimensions of operators governing the critical fluctuations (related to critical 
exponents) are evaluated from the finite- A" corrections of the excited states. 
For each primary operator, with dimension x^, in the operator algebra of the 
system, there exists an infinite tower of eigenstates of H whose energy m , 
are given by 

= ei)1 + _£(^ + m + m') + o(iV- 1 ), (7) 

where to, to' = 0, 1, . . . . The real part of the energies appears in (|7|) since 
in general they are complex due to the non-hermicity of H. The multiplicity 
of energies with the same real part 3? (e^ m /(A^)^ are given in terms of the 
product of two Virasoro algebras, as in the periodic case. We tested exten- 
sively the relations (§) and (0) for several excited eigenenergies in the case 
of the Ising model with helical boundary conditions. 



9fc (c 
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3 The spin-1 Model 



In order to test our numerical methods and the conformal invariant relations 
(ID and (0), we consider initially the Blume-Capel model with spin S = 1. 
Earlier finite-size studies that explores the conformal invariance of this model 
were done in its r-continuum formulation || pLOfl . In figure 1 we show the 
extrapolated (N — > oo) transition curve t c = t c (d) (continuum line in the 
figure) obtained by solving Eq. for lattice sizes up to iV = 12. Note that 
for d <C — 1 only configurations where Sj = ±1 is allowed and we have an 
effective Ising model with t c = 2/ ln(>/2 + 1) = 2.269.... 

The expected tricritical point of the model is obtained by extrapolating 
the sequences obtained by solving Eq. (|J). Our estimate for this point is 
t t = 0.609(4) and d t = 1.965(5) (we solved (D for lattice sizes up to iV = 10), 
in agreement with previous Monte Carlo and finite-size calculations ||. 
In [U these estimates were done from data of lattice sizes up to iV = 8. 
The extrapolations along this paper were obtained by using the e-alternated 
VBS approximation [p4}| . The errors are estimated from the stability of the 
extrapolations and are in the last digit. In figure 1 we show in a large scale 
the region where the tricritical point is located. In this figure this point is 
the end point of the continuum curve. 

Strictly speaking the relation (f|) only gives continuous phase transitions. 
At the first-order line, in the spin-1 model, we have the coexistence of three 
phases: two ordered ferromagnetic (< Si >^ 0) and one disordered phase 
(< S{ >= 0). Consequently the gaps related with the three lowest eigenval- 
ues vanishes exponentially as the lattice size increases. A possible finite-size 
estimate for the first-order transition line, which we test in this paper, is ob- 
tained by the following procedure. For a given lattice size N we calculate the 
points where the gap corresponding to the third eigenvalue has its minimum 
value. The extrapolation iV —>■ oo of these points give us our estimate for 
the first-order transition line. The dotted line shown in Fig. 1 was obtained 
by this procedure. As we can see in this figure, the estimated first-order 
transition line (dotted line) finishes at the tricritical point. 

In the critical regions of the phase transition line (continuum curve) the 
conformal anomaly and the scaling dimensions can be calculated exploring 
the conformal invariant relations (|6|) and (|7]). From Eq. (|6|) a possible way 
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to extract c is by extrapolating the sequence 



r N,N + i _ 6 ( e lt i(N+l) _ ei,i(iV) \ / J_ _J__\ ^ 

"tt^ iV+1 N )\Ni (N + iy) ' lJ 

calculated at t c (d). Examples of such sequences together with the extrapo- 
lated results for the spin-1 model, are shown in the first four columns of table 
1. We clearly see from this table the expected results, i. e., an Ising behavior 
with conformal anomaly c = 1/2 in the left of tricritical point (continuum 
curve) and c = 7/10 at the tricritical point. 

From Eq. (0) the scaling dimensions x(n, p) related to the nth (n = 1,2...) 
energy in the sector with parity p can be obtained by extrapolating the 
sequence 

xN = WvWl-uW)^ (9) 

2n 

The smallest dimension x(l, —1) corresponds to a Z(2)-order parameter. 
Some of the related finite-size sequences are presented in the first four columns 
of the table 2 for the spin-1 model. As expect this dimension has the value 
1/8 for the critical points in the Ising universality class and 3/40=0.075 at 
the tricritical point. The conformal towers at tricritical point was obtained 
previously by using the r-continuum Hamiltonian formulation of the model 

In order to compare our estimate for the tricritical point with the previous 
results we have also calculated the conformal anomaly and scaling dimensions 
by using the tricritical point estimated in Refs. [0, |8|. These results show 
us that the extrapolated values of the conformal anomaly is similar to ours, 
however our estimate for the scaling dimensions are much better than the 
values obtained at their estimated tricritical point. For these reasons we 
believe that the tricritical point we found is a better estimate. 



4 The spin-3/2 Model 

We now apply the numerical methods of section 2 in the controversial case 
of spin 3/2. In figure 2 we show the extrapolated transition curve t c = t c (d) 
(continuum line), obtained by solving Eq. (f|) for lattice sizes up to iV = 10. 
The limiting values t c = 9/(21n( v / 2 + 1) = 5.106... (d < -1) and t c = 
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l/(21n(v / 2 + 1) = 0. 567.. .(d ^> 1) can be easily understood since in these 
limits the model reduces to an Ising model. 

As we mentioned in the introduction there is a controversy in the lit- 
erature concerning the existence or not of a multicritical point along the 
transition line. In order to clarify this point we try to obtain finite-size es- 
timates for this point by using relation (|5|). Our results however show no 
consistent solutions of relations @ for lattice sizes up to iV = 10, which in- 
dicate the absence of a multicritical point along the transition line. Another 



method, also used to locate multicritical points p2| , is obtained from the si- 



multaneous crossing of two different gaps on a given pair of lattices (instead 
of three lattices as in (f|)). Trying several different gaps we also did not find, 
within this method, a multicritical point along the transition curve. 

The first-order transition line was also estimated in a similar way as we did 
for the spin-1 model in last section. Along this line we have the coexistence 
of four ordered phases, with most of spins in state s = ±3/2, 1/2 or —1/2 
respectively. The finite-size estimator, for a given N, is given by the points 
where the fourth gap has its minimum value, since in the bulk limit the four 
lowest eigenenergies will degenerate. The extrapolated curve (N — > oo) is 
the dotted curve shown in Fig. 2. As depicted in a large scale in this figure, 
this first-order transition line does not touch the continuous transition line 
(differently from the case S — 1), in agreement with the topology predicted 
for the phase diagram by mean field calculations [15| . 



In order to supplement our results by an independent method we also did 
Monte Carlo simulations for the spin-3/2 model. Details of the simulations 
are presented in appendix A. The points for the phase diagram are the squares 
shown in Fig. 2. Taking into account the scale of this figure the agreement of 
the Monte Carlo simulations (squares) and the transfer matrix results (circle) 
is good. These results indicate the absence of a multicritical point along the 
continuous transition line. According to this scenario the whole continuous 
phase-transition line should belong to the Ising universality class. In order 
to illustrate this point we present in the last four columns of table 1, for 
some values of d, the finite-size sequences (|8|) of the conformal anomaly of 
the model. In the last four columns of table 2 we also show some of finite-size 
sequences (§) for the dimension x(l, —1) of the order parameter. Again we 
obtain the expected results 1/8 of the Ising-order parameter. 

As usual in the region where the first-order phase transition takes place 
the finite-size sequences have slower convergence due to crossover effects. 
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This fact prevent us to get a good estimate for the endpoint of the first 
order transition line. Our Monte Carlo simulations ( see appendix A) give 
an estimate t I = 0.80 and dj = 1.96, for this point. 



5 The r-continuum quantum Hamiltonian of 
the spin-3/2 model and other related quan- 
tum chains 

Our results in the last section indicate that the spin-3/2 Blume-Capel model 
does not have a multicritical point along the transition line that separates 
the disordered phase from the ordered ones. Instead of working in the Eu- 
clidean version, we can also consider the model in its r-continuum Hamilto- 
nian formulation. This quantum Hamiltonian which is obtained by imposing 
extreme anisotropic relations between its couplings in the horizontal and ver- 
tical directions, can be derived straightforwardly (see jn|) for example). For 
periodic boundary conditions it is given by H(b — 1) where 



N 



H{b) = ~ E SISl+i - MS?) 2 ~ x 2 R i (b) 1 



(10) 



and Xi and X2 are coupling constants. The operators Si and Ri(b) in the 
basis where Si is diagonal are given by 

Sf = 1 <g> 1 <g> • • • 1 g) S z ® 1 • • • <g> 1, 

R,(b) = 1 ® 1 ® • • - 1 ® R(b) ® 1 • • • ® 1, 

where the matrices S z and R(b) are in the ith position in the product and 
are given by 



R 
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Diagonalizing the Hamiltonian ( |10D for several lattices we try to locate a 
multicritical point, by solving Eq. (f|). In agreement with the absence of 
multicritical behavior, we found no triple crossing for lattice sizes N > A. 
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It is important to mention that, on the contrary to these results, the 
spin- 3/2 Hamiltonian H(2/y/3), given by (|10D is known to have a tetracritical 
point along the transition line separating the disordered and ordered line [E3 



|25fl . The matrices S z and R(2/y/3), given in (|IID , are in this case the standard 
^-diagonal representation of the spin-3/2 SU(2) matrices. This tetracritical 
point was located |23| by the same techniques we used in this paper and 
its long-distance physics is related by a conformal field theory with central 
charge c = 4/5. Its is interesting to observe that such an small difference 
in the Hamiltonians produce such differences in the critical behavior of the 
Hamiltonians. In order to better see this fact numerically, we study the 
Hamiltonian ( |10D with b ranging from 1 to 2/y / 3- We verify that solutions 
of Eq. (7), for multicritical points happens only for values of b bigger than 
or very close to 2/y/3. 

These results in the quantum Hamiltonian formulation indicate that mul- 
ticritical points along the order-disorder line can be generated only if we in- 
troduce additional interactions in the spin-3/2 Blume-Capel model. The net 
effect of the coupling b in the quantum chain (|T0D is to control the preference 
for neighboring spins (s, s') along the time (vertical) direction to be in pairs 
(1/2,-1/2) or (-1/2,1/2). One of possible ways to introduce this effect in the 
classical model is to add in the spin-3/2 Blume-Capel ([!]) Hamiltonian an 
extra term V s>s '(a) = — a(5 S) i/ 2 5 s ',--i/2 + ^',1/2^,-1/2), along all the links in 
the vertical. The product (3a will play a similar role as b in ([L0|). Studying 
numerically the transfer matrix of such classical Hamiltonian we found for 
a = 1 a multicritical point at t t = 0.87(8) and d t = 1.60(5). Our calcula- 
tion at this point are in favor of a value c = 4/5 and x(l, —1) = 0.05 for 
the conformal anomaly and the dimension of most relevant order parameter, 
respectively. These are the predicted values for a minimal conformal theory 
with m = 5 p5| . 



6 Summary and Conclusion 

Our aim in this paper was the study of the critical properties of the spin-3/2 
Blume-Capel model by using the methods of finite-size scaling in strip ge- 
ometries, conformal invariance and Monte Carlo simulations. The methods 
and useful relations for these calculation were presented in section 2 and in 
appendix A, and were tested for spin-1 Blume-Capel Model with a precise 
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evaluation of its phase diagram. Our results for the spin-3/2 help to choose 
between the two existing contradictory predictions for its phase diagram in 
(t, (i)-parameter space. Mean-field calculations |T5| predict (scenario a) a 



phase diagram where there is no multicritical point along the second-order 
phase transition line. In this scenario a multicritical point will exist as an 
isolated end point of a first-order transition line. On the other hand approx- 
imate calculations based on renormalization- group methods in real space 



11], |12| give us (scenario b) a multicritical point along the disorder-order 
transition line. This point being the end point of the first-order transition 
line. Since the mean-field approximation is expected to be correct only for 
high dimensions and the other calculations predicting scenario b were done 
directly in two dimensions we would expect, at least for two dimensions sce- 
nario b as the correct one. In parallel to these results there exist in the 
literature J23, precise calculations on a related spin-3/2 quantum chain, 
which also indicate a multicritical point along the transition line. Although 
these calculations were done in a quantum Hamiltonian these results induce 
us in favor of scenario 6, since this quantum version chain (H(2/y/3) in Eq. 
(|10D ) has naively the same type of couplings as the spin-3/2 Blume-Capel 
Model. 

However, to our surprise, our results presented in section 4 indicate the 
scenario a as the correct one for the spin-3/2 Blume-Capel Model. In section 
5 we also analyze this last model in its r-continuum Hamiltonian formula- 
tion {H{1) in Eq. (|i~0|) ) for the sake of comparison with the above mentioned 
spin-3/2 quantum chain (H(2/\^3) in Eq. (0)), which shows multicritical be- 
havior along the order- disorder transition line. Although these Hamiltonians 
are very similar our finite-size scaling studies indicate that these differences 
are enough to produce quite distinct critical behavior. The results of section 
5 also indicate that a phase diagram for a spin-3/2 model like that in scenario 
b can be produced by adding extra terms in the Blume-Capel model. 
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Appendix 



A Monte Carlo Simulations 

Our Monte Carlo simulations for the spin-3/2 Blume-Capel model were done 
on square lattices with N 2 sites and periodic boundary conditions in all 
directions. We have used the standard single spin flip algorithm [EB] where 
the transition probability from a system configuration with energy Ei 
(where a chosen spin is in a state i) to another configuration with energy Ej 
(where this spin is now in a state j) is 

where AE^ = Ej — E^. Since the spin variable can assume more than two 
states we have chosen, at random, the state j from the (2S — 1) possible 
distinct states from state i. 

The thermodynamic quantities computed in these simulations are the 
magnetization per site 

m = (M) , 

the magnetic susceptibility 

Xm = (3N 2 ((M-m) 2 } , 

where 

M EiSl 



N 2 ' 

and the quadrupole magnetic susceptibility 

Xq = PN 2 ((q-(q)) 2 ) , 

where 

Ys 2 
q N 2 

Typical runnings have been made on a periodic lattice with N = 64 spins 
and included 6.0 x 10 4 iterations or Monte Carlo Steps (MCS) per spin of 
the algorithm and the averages of the measured quantities were taken after 
discarding the first 1.0 x 10 4 MCS per spin. We have found that by changing 
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the lattice size from N = 64 to iV = 128 the relevant measured quantities did 
not change appreciably. For this reason, in most of the computations we have 
taken N = 64 in order to save computing time. In addition, we have also 
calculated some Monte Carlo averages taking data from: i) every MCS per 
spin; ii) every 10 MCS per spin; and iii) runnings including 1.2 x 10 5 MCS per 
spin with averages taken after discarding the first 2.0 x 10 4 MCS per spin. 
In all these runnings, no significant differences have also been observed in 
the thermodynamic quantities of interest when compared to the typical one 
previously stated in the beginning of this paragraph. These results assure 
us that possible correlations among different configurations, if any, are not 
relevant for the computation of the magnetization and susceptibilities. In- 
deed, our results, within the present numerical precision, are good enough to 
be compared to other values obtained from different methods. For instance, 
at d = we have t = 3.32 and in the spin- 1/2 Ising limit d — > oo we have 
t = 0.58 which is comparable to the exact value t = 0.57. 

The simulations for the phase diagram close to the isolated multiphase 
critical point are shown in Fig. |2|. Second-order phase transitions have been 
determined by the strong peak in the magnetic susceptibility Xm, while first- 
order transitions have been located through the discontinuous behavior of the 
magnetization accompanied by a strong peak in the quadrupole susceptibility 
X q - A typical example is shown in Figs. [3], |4] for d = 1.998 as a function of the 
reduced temperature t. Such behavior of Xq is useful for detecting first-order 
transitions when the discontinuity in the magnetization of the two ordered 
phases is very small. The discrepancy between the Monte Carlo simulations 
and finite-size calculations for d very close to 2 reflects the large metastability 
of the ordered phases. For instance, a metastable phase with m = 3/2 will 
persist at low temperatures for d slightly above 2 due to the ordering state 
starting configuration (m = 3/2). In fact, by taking a random-state starting 
configuration and sweeping up and down the temperature for constant d > 2 
the only stable phase found is the one having m— 1/2 for T = 0. The point 
at d = 2.0 and t = 0.4 has not been simulated and is in Fig. 2 just for the 
sake of clarity. 

The isolated multiphase critical point in the ordered phase is rather diffi- 
cult to be located with reasonable precision through Monte Carlo simulations 
(an estimate can be given by dj = 1.96 and tj = 0.80). Despite of that, we 
have carefully swept paths of constant temperatures for t > 0.81 and constant 
d for d > 1.92 and we have not found any indication of the first-order line, 
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separating the low temperature ordered phases, terminating in the second 
order transition line. This is in disagreement with previous renormalization 
group calculations [O, ITS] and reproduce the qualitative behavior predicted 



by mean field like calculations [15], and are in agreement with our results of 



finite-sizes calculations on strips (see Fig. 2). 
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Figure and Table Captions 



Figure 1: Phase diagram in the (d — t)-plane for the spin-1 Blume-Capel 
model. The continuum line is the second-order transition line and the dotted 
line is the first-order transition. The tricritical point is located at the end 
point of the continuum line. 



Figure 2: Phase diagram in the (d — t)-plane for the spin-3/2 Blume-Capel 
model. The continuum line is the second-order transition line and the dotted 
line is the first-order transition. The circles were obtained by finite-size 
estimates and the squares by Monte Carlo. The full square represents the 
corresponding isolated multiphase critical point. There are two different 
ferromagnetic phases with mi — > 3/2, and rri2 ^l/2asT-^0. The open 
squares in the blown-up region 1.90 < d < 2.10 are obtained by Monte Carlo 
simulations. 



Figure 3: Quadrupolar susceptibility (open squares) and magnetization 
(open circles) as a function of t for the Blume-Capel model with S = 3/2 and 
d = 1.998. Xq nas a strong peak only close to t — 0.52 and the magnetization 
a discontinuity of about 0.94, indicating a first-order transition. The data 
were taken for increasing t. The lines are guide to the eyes. The data 
have been normalized: the maximum magnetization is m = 1.499 and the 
maximum susceptibility is Xq — 2.02. 
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Figure 4: Magnetic susceptibility (open squares) and magnetization (open 
circles) as a function of t for the Blume-Capel model with 5 = 3/2 and 
d = 1.998. Xm has a strong peak only close to t — 0.66 and the magnetization 
is continuous, indicating a second-order transition. The data were taken for 
increasing t. The lines are guide to the eyes. The data have been normalized: 
the maximum magnetization is m = 1.499 and the maximum susceptibility 
is Xm = 53.6. 



Table 1 - Finite-size data c N,N+l given by (§) and extrapolations for the 
conformal anomaly for some values of d. The first four columns are the 
values for the S = 1 model while the last four are those of the S = 3/2 
model. 



Table 2 - Finite-size data x N (l, —1) given by @ and extrapolations for the 
scaling dimension of the parameter order for some values of d. The first four 
columns are the values for the S = 1 model while the last four are those of 
the S = 3/2 model. 
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Table 1 



N 


S=l 


S=3/2 


d=-0.5 


(2=0 


d=1.9 


ci=1.9655 


d=-0.5 


(2=0 


c2=0.5 


(2=2 


t c =l. 794(7) 


t c =1.681(5) 


t c =0.764(7) 


2 C =0.609(4) 


t c =3. 549(8) 


t c =3.287(2) 


t c =2. 972(3) 


2 C =0.64(7) 


4 


0.460931 


0.461723 


0.570695 


0.641586 


0.461970 


0.463160 


0.465282 


0.626807 


5 


0.478008 


0.478429 


0.578328 


0.666641 


0.478574 


0.479283 


0.480650 


0.611765 


6 


0.486062 


0.486307 


0.575809 


0.678638 


0.486400 


0.486858 


0.487798 


0.589182 


7 


0.490323 


0.490480 


0.569959 


0.685059 


0.490544 


0.490864 


0.491549 


0.568249 


8 


0.492838 


0.492948 


0.563332 


0.688878 


0.492994 


0.493232 


0.493754 


0.551328 


9 


0.494458 


0.494539 


0.556880 


0.691349 


0.494573 


0.494758 


0.495170 


0.538405 


10 


0.495569 


0.495633 


0.550943 


0.693048 


0.495659 


0.495807 


0.496142 


0.528798 


11 


0.496370 


0.496421 


0.545616 


0.694269 


0.496440 


0.496562 


0.496841 


0.521749 


12 


0.496968 


0.497009 


0.540896 


0.695175 










13 


0.497427 


0.497462 


0.536738 


0.695864 


14 


0.497788 


0.497817 


0.533085 


0.696401 


oo 


0.499(9) 


0.499(9) 


0.50(3) 


0.698(9) 


0.499(8) 


0.499(8) 


0.499(8) 


0.50(5) 



Table 2 



N\ d 


S=l 


S=3/2 


-0.5 





1.9 


1.9655 


-0.5 





0.5 


2 


4 


0.118424 


0.118238 


0.100130 


0.071623 


0.118424 


0.117745 


0.117610 


0.089071 


5 


0.121069 


0.120923 


0.105237 


0.073095 


0.121185 


0.120449 


0.120484 


0.100904 


6 


0.122389 


0.122278 


0.108791 


0.073783 


0.122612 


0.121806 


0.121985 


0.108882 


7 


0.123137 


0.123049 


0.111422 


0.074134 


0.123454 


0.122565 


0.122863 


0.113933 


8 


0.123602 


0.123532 


0.113439 


0.074315 


0.124003 


0.123022 


0.123424 


0.117077 


9 


0.123913 


0.123855 


0.115022 


0.074400 


0.124391 


0.123312 


0.123807 


0.119047 


10 


0.124132 


0.124082 


0.116285 


0.074425 


0.124683 


0.123500 


0.124082 


0.120304 


11 


0.124292 


0.124248 


0.117305 


0.074409 


0.124913 


0.123625 


0.124287 


0.121121 


12 


0.124413 


0.124373 


0.118136 


0.074362 


0.125102 


0.123706 


0.124446 


0.121662 


13 


0.124507 


0.124470 


0.118821 


0.074290 










14 


0.124582 


0.124547 


0.119387 


0.074199 


15 


0.124642 


0.124608 


0.119860 


0.074091 


oo 


0.1250(0) 


0.1250(2) 


0.12(2) 


0.074(0) 


0.13(0) 


0.123(8) 


0.125(5) 


0.122(2) 
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